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Abstract 

Two possible diagnostics of stretching and folding (S&F) in fluid flows 
are discussed, based on the dynamics of the gradient of potential vor- 
ticity (q = u> ■ V6) associated with solutions of the three-dimensional 
Euler and Navier-Stokes equations. The vector B = Vg x V0 satisfies 
the same type of stretching and folding equation as that for the vorticity 
field us in the incompressible Euler equations (Gibbon & Holm, 2010). 
The quantity 9 may be chosen as the potential temperature for the strat- 
ified, rotating Euler/Navier-Stokes equations, or it may play the role of 
a seeded passive scalar for the Euler equations alone. The first discus- 
sion of these S&F-flow diagnostics concerns a numerical test for Euler 
codes and also includes a connection with the two-dimensional surface 
quasi-geostrophic equations. The second S&F-flow diagnostic concerns 
the evolution of the Lamb vector D = u> x u, which is the nonlinearity 
for Euler's equations apart from the pressure. The curl of the Lamb 
vector (zu := curlD) turns out to possess similar stretching and folding 
properties to that of the B-vector. 



PACS numbers: 47.10.A-, 47.15.ki 



1.1 Introduction 

This paper considers two variants of the stretching and folding proper- 
ties of gradients of solutions of the three-dimensional Euler and Navier- 
Stokcs equations, following recent work of Gibbon & Holm (2010). Fine- 
scale structures, diagnosed in either inviscid or viscous turbulence and 
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MHD by the presence of large gradients, are created in the tortuous 
stretching and folding processes that arise from the vortex stretching 
term in the Euler and Navier-Stokes equations. These fine-scale struc- 
tures are not wholly understood, as they lie at the heart of unsolved 
regularity issues that have challenged mathematicians for more than a 
generation. 

In what follows the advected scalar field 9 will comprise either : (i) the 
potential temperature for the stratified, rotating Euler/Navier-Stokes 
equations ; or (ii) a passive scalar for the Euler equations alone. The 
main theme revolves around the role of the vector 3 = Vq x V6* where 
the potential vorticity q = u> ■ V9 is conserved on fluid particle paths 
in either case. The basis of the result, already discussed by Kurgansky 
& Tatarskaya (1987) and Kurgansky & Pisnichcnko (2000), is that in 
the incompressible Euler case the vector 3 satisfies the same equation 
as that for the vorticity, thus suggesting intriguing stretching and fold- 
ing properties for the gradient of the projection of u) on the normal to 
level surfaces of 9. This result, summarized in HI. 21 also has interesting 
consequences for the Navier-Stokes equations (Gibbon & Holm, 2010). 

The first of the variants on this theme in HI. 3. II concerns a scheme for 
testing numerical Euler codes which have been designed to address the 
issue of whether the equations develop a finite time singularity. Thus 
it is apposite to devote the introduction to this section HI. 31 to listing 
some of the Euler literature in this area. A closely associated prob- 
lem forms the subject of HI. 3. 21 in which a connection is established 
with the two-dimensional surface quasi-geostrophic equations (2D-QG) 
studied by Constantin, Majda & Tabak (1994). In this X 7± 9 in two- 
dimensions obeys the same vortex stretching equation as that of uj for 
three-dimensional Euler. It is shown that the 2D-QG equations are em- 
bedded as a special case in the equation for 3. 

The second main variant revolves around the stretching and folding 
properties of the Lamb vector D = u> x u for the incompressible Euler 
equations. The Lamb vector comprises the nonlincarity of the Euler 
equations aside from the pressure, so its evolution is of importance. In 
§?? it is shown that its curl 

zu := curl.D = curl (uj x u) (1-1) 

also satisfies the same type of stretching equation as that for 3, while 
its divergence (divD) obeys a continuity equation. In the compressible 
case this may have interesting consequences for the study of jet-noise 
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although this is beyond the present scope of this paper. Section 11.4.21 
deals with the evolution of the gradient of helicity density A = w • u 
which also appears to possess similar stretching and folding properties. 

Let us begin with the notation for the incompressible Euler equations, 
which arc expressed as 

Du 

— = -Vp, divu = 0, (1.2) 



d t u-uxuj = -V(p+±u 2 ). (1.3) 

The chosen domain is a three-dimensional periodic box V = [0,L] 3 . u 
is the velocity field of the fluid and the material derivative is defined by 

^- t =d t +u-V. (1.4) 
The vorticity field us = curl u satisfies 

dtU) — curl(tt xw) = 0. (1.5) 
This formula can also be written in the familiar vortex stretching format 

— = W'VmeSw, (1.6) 

where Sy = \ ( U U + u j-0 ^ s the rate of strain matrix. Equations (11.51) 
and (|1.6[) are equivalent evolution equations for oj. Euler data roughens 
very quickly, a fact which is mainly due to the stretching and folding 
processes caused by the rapid alignment or anti-alignment of u> with 
positive and negative eigenvectors of S. 

The main aim of this paper is to show that these stretching and folding 
processes are shared by several other variables in the Euler and Navier- 
Stokcs equations. 

While the existence of some very weak solutions has been proved 
(Shnirelman, 1997; Brenier, 1999; Majda & Bertozzi, 2001; De Lellis 
& Szekelyhidi, 2007, 2008; Brenier, De Lellis & Szekelyhidi, 2009), nev- 
ertheless Lcray-type weak solutions arc unknown. However, if we are to 
progress in our understanding of the properties of solutions of the Euler 
equations, our lack of knowledge forces us to make some assumptions 
about the existence of solutions. The fundamental result on existence 
of solutions of the three-dimensional Euler equations is the theorem due 
to Beale, Kato & Majda (1984) which is assumed to hold : 
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Theorem 1.1.1 (Beale, Koto & Majda, 1984) There exists a global 
solution u <G C([0, oo];H s ) D C 1 ([0, oo];i? s_1 ) of the Euler equations 
for s > 3 if and only if, for every t > 0, 

t 

rJHoodr < oo. (1.7) 

Remarks, (i) The value of this result is that computationally only the 
quantity ||w||oo needs to be monitored. If this is finite everywhere in the 
domain of flow at a time t then the solutions are regular at that time. 

(ii) It does not predict a singularity in ||w||oo but it restricts those that 
may potentially occur of the type ||tt>||oo ~ (t s — t)~ p to the range p > 1. 
When p < 1 the theorem is violated. 

(iii) Kozono and Taniuchi (2000) have proved a version of this theorem 
in the BMO-norm (bounded mean oscillations) which is slightly weaker 
than the L°°-norm. 

(iv) Further analytical approaches have centred around conditional esti- 
mates on the magnitude and direction of vorticity that include the direc- 
tion of vorticity. These are extensions of the Beale-Kato-Majda theorem ; 
the most significant papers arc those by Constantin, Fcfferman & Majda 
(1996), Deng, Hou & Yu (2005,2006) and Chae (2003,2004,2005,2007). 



1.2 The 3 equation for the stratified Euler &: Navier-Stokes 

equations 

1.2.1 The stratified, rotating Euler equations 

Let us consider the three-dimensional incompressible Euler equations 
for an incompressible, stratified, rotating flow (f2 = k f2) in terms of the 
velocity field u(x, t) and the potential temperature 9 

Du 

— + 2(nxu) + a k6 = -Vp, (1.8) 

where ao is a dimcnsionless constant and where 0(x, t) evolves passively 
according to 

^=0. (1.9) 

Dt y ' 

How 8(x, t) and other variables might accumulate into large local con- 
centrations is of intercst0. To pursue this, consider the vorticity u> = 

1 The BKM theorem expressed in the last section is valid when 8 is no more than 
a passive scalar driven by an Euler flow as in 11.51 . For stratified Euler 1 11 .86 together 
with 11.91 . however, it is necessary to assume that f£ (||a;||oo + ||V0||oo) d-r is finite. 
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curlu and define u> ro t = oj + 2J~2, which satisfies 

= u, rot ■ Vu - V^9 V 1 - = (d y , -d x , 0) . (1.10) 

The potential vorticity defined by (Hoskins, Mclntyre & Robertson, 
1985) 

q = oj rot ■ V6, (1.11) 

satisfies Ertcl's theorem (Ertcl, 1942; Truesdell & Toupin, 1960; Ohki- 
tani, 1993; Kuznctsov & Zakharov, 1997) because 

U) rot ■ VU ■ VV + UJrot ■ V 



Dt V D t J \Dt J 

= -\/ ± e-ve = o. (i.i2) 




= const 



Fig. 1.1. The vector B points along an intersection of level sets of the two 
Lagrangian flow constants (q, 9) 



This establishes two quantities q and 9 that are each conserved along 
flow lines, and whose level sets intersect as in Figure 11.11 Then with B 
defined as 

B = \7qxV6 (1.13) 

Kurgansky & Tartskaya (1987), Kurgansky & Pisnichenko (2000) have 
observed that B satisfied 

d t B = curl (uxB) . (1.14) 

x Notc that this is a Clebsch representation of the divergence- free vector 0, not a 
decomposition of the vorticity oj. See Ohkitani (2008) for a recent study of the latter 
and Holm & Kupershmidt (1983) for a review of the Clebsch approach. Moreover, 
the helicity of B given by f n B ■ curl _1 jBdy is necessarily zero for homogeneous or 
periodic boundary conditions. 
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Of course this may be written equivalently in the familiar vortex stretch- 
ing format f| 1 .€>[) 

^ = B-V« (1.15) 

thereby highlighting the fact that alignment of 3 with eigenvectors of 
Vit is critical to the stretching process. In Figure 11.11 the vector 3 
is tangent to the curve defined by the intersection of q = const and 
9 = const. Thus, 3 plays the same role as that for iv and for the 
magnetic B-field in MHD (Moffatt, 1978; Palmer, 1988). Hence, all the 
stretching and folding properties associated with vorticity or magnetic 
field-lines also apply to 3 even though 3 contains various projections of 
u>, Vu, V8 and VV#. Moreover, for any surface S(u) moving with the 
flow m, one finds 

^ [ 3-dS = 0. (1.16) 

S(u) 



dt 



1.2.2 The stratified Navier- Stokes equations 

Now let us turn to the Navier-Stokes equations coupled to the 0-ficld. 
(In what follows the rotation will be ignored.) These equations are 

Du 

— + a 6k = Re^Au-Vp, (1.17) 

^ = (aRe)- 1 A6. (1.18) 

Here, the potential vorticity q = ui ■ V9 is no longer a material constant 
but, instead, evolves according to 

Dq ( D " „ \ ™ ^ (DO 

— = w • Vu • V(9 + oj • V — 

Dt \Dt J \Dt 

= (Re- 1 ALJ-V ± 6) ■V6' + a;-V[(cri?e)- 1 A6l] 

= div{iie _1 Att X V0 + (aRe^uiAe} . (1.19) 

The material advection property is destroyed but the introduction of a 
transport velocity field U q transforms (|1.19|) into a continuity equation 

d t q + div (qU q ) = , (1.20) 

thus making q a PV density, and where tl q is defined through the relation 

q(U q -u) = -Rer 1 (Au x V6> + a^uAO) . (1.21) 
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The introduction of the transport velocity held U q is originally due to 
Haynes & Mclntyre (1987). Note that d\vU q ^ although d\vU q = 
O (Re) 1 . Consistent with numerical studies on reconnection (Herring, 
Kerr & Rotunno 1994), this scaling with Reynolds number Re may sug- 
gest that in the early or intermediate stages of a flow this divergence 
may be small. What about the evolution of the variable 0? It is easily 
seen that 

d t Q + U q -VQ = d t 6 + u- VO-Re^q- 1 {Au x V0 + er _x wA0} • V0 
= d t 8 + u- V0- (aRe)' 1 AO = . (1.22) 
The formal result for the stratified Navier-Stokes equation is : 

Theorem 1.2.1 The scalar quantities q and satisfy 

d t q+ div(qU q ) = 0, d t e + U q -V6 = 0, (1.23) 

and 3 = Vq x V0 satisfies the stretching and folding relation 

d t B - curl(U q x S) = V q , (1.24) 

where the divergence-less vector T) q is given by 

V q = -V(q divU q ) x V0 , (1.25) 

and the transport velocity tl q is defined as in (1.21]) . Moreover, for any 
surface S(M q ) moving with the flow U q 

f B dS= [ T> q -dS. (1.26) 

dt JS(u q ) JS(u q ) 

This is the natural way of expressing problems in the vortex stretching 
format using the transport velocity U q . 

1.3 The Euler singularity problem 

Out of large-scale computations of solutions of the three dimensional Eu- 
ler equations has emerged the natural question of whether a singularity 
develops in a finite time (Majda & Bertozzi, 2001; Bardos & Titi, 2007; 
Constantin, 2008; Gibbon, 2008). An extensive literature has arisen on 
this question but no conclusion has yet been agreed : sec Bardos & Bcna- 
chour (1977); Morf, Orszag & Frisch (1980); Chorin (1982), Brachct, Me- 
iron, Orszag, Nickel, Morf & Frisch (1983); Siggia (1984); Kida (1985); 
Ashurst & Meiron (1987); Pumir & Kerr, (1987); Pumir & Siggia (1990); 
Grauer & Sidcris (1991); Bell & Marcus (1992); Brachct, Mcneguzzi, 
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Vincent, Politano & Sulem (1992); Kerr (1993, 2005a, 2005b); Boratav 
& Pelz (1994,1995); Pelz (1997,2001); Pelz & Gulak (1997); Grauer, 
Marliani & Germaschewski (1998); Cichowlas & Brachet (2005); Gu- 
lak & Pelz (2005); Pelz & Ohkitani (2005); Pauls, Matsumoto, Frisch & 
Bee (2006). Regarding more recent work, the fine-scale computations by 
Hou and Li (2006, 2007) that see only super-exponential growth in us con- 
tradict both the older computations of Kerr (1993, 2005a) together with 
newer results by Bustamante and Kerr (2007), in which a finite time sin- 
gularity has been observed. Similar but not identical anti-parallel vortex 
tube initial conditions have been used in the two bodies of results which 
largely coincide until a late stage. Two further recent contributions are 
those of Orlandi and Carnevale (2007) who used initial conditions in 
the form of Lamb dipoles to observe singular behaviour, as did Grafke, 
Homann, Dreher & Grauer (2007). 



1.3.1 A numerical test for Euler computations 

Let us take the Euler equations for the velocity field it in their standard 
form as in (|1.2[) without the rotation or buoyancy used in m.2\ The new 
feature of the proposed test is to introduce a passive tracer concentration 
6(x, t) satisfying 

and whose initial data arc under the investigator's control. Introducing 
allows us, as before, to use q = uj ■ V6> which still obeys 

§f=0, (1.28) 

and which therefore allows the use of the same definition B = Vq x V0. 
This is endowed with initial conditions inherited from those for it and 
9. B must satisfy 

DB 

d t B = curl (it x B) or — = B ■ Vit , (1.29) 

which mimics the vorticity stretching equation (|1.5[) . The critical point 
about the use of the vector field B is that it has embedded information 
on u>, Vu;, V# and VV#. It evolves in the same way as a> and so is 
subjected to similar stretching and folding processes. It can, however, 
be evaluated at any particular time t in several distinct ways : it can 
be evaluated from the result of the evolution in (|1.29[) at time t, or it 
can be computed from its definition using it and 6 evolved up to and 
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evaluated at time t. The degree to which these distinct evaluations 
agree or disagree provides a quantitative gauge of the accuracy of the 
numerical computation. It is not clear that there is a natural scale for 
the inevitable discrepancies produced in any particular computation. 
However, this procedure produces a precise diagnostic quantity that, 
given identical initial data, can be directly compared side-by-side for 
different numerical computations to evaluate their relative accuracy. The 
suggested test is : 

(i) Choose initial data for u and 9, thereby fixing initial data for q 
and B. 

(ii) Evolve u and simultaneously solve Dd/Dt = 0, Dq/Dt = and 
DB/Dt = B • Vm. 

(hi) Test the resolution at any time t > by constructing <ji(- , t) = 
u){- , t) ■ V0(- , t) and then: 

(a) compare the solution for B(-, t) obtained from solving the 
stretching equation DB/Dt = B ■ Vm with 

8i(-,t) = Vgi(-,t)xV0(-,f) (1.30) 

(b) Furthermore, compare this with 

B 2 (- , t) = Vg (• , t) x V0(- , t) (1.31) 

where q (• , t) is the evolved solution of Dq/Dt = 0. 

(iv) For fixed initial data for u this procedure may be implemented for 
a variety of "markers" 9 n (- , t) evolving from distinct initial data 
9 n (- , 0) to diagnose the numerical accuracy in different regions 
of the flow. 

Because B contains Vu, comparing the different computations of B, B\ 
and 2?2 tests the accuracy of the computation of some of the small scale 
structures in the flow. Given the generally acknowledged difficulties in 
accurately computing the evolution of passive scalars such as in (|1.27l) 
or (|1.28p , how initial data are chosen may be critical to the calculation. 
In particular, the appearance of null points in the vorticity field may 
create significant obstacles - see Ohkitani (2008). 

1.3.2 Connection with the two-dimensional quasi- geostrophic 

equations 

The open nature of whether the three dimensional Euler equations de- 
velops a singularity naturally leads to the idea of studying other simpler 
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problems that mimic this behaviour. The foremost example is the case of 
the two dimensional surface quasi-geostrophic (2D-QG) equations. The 
strong fronts observed in numerical computations though the existence 
of a vortex stretching term have led Constantin, Majda & Tabak (1994) 
to suggest that these might model singularity development in the three 
dimensional Euler equations. It turns out that the 2D-QG equations are 
embedded in the equation for B in the following way. 

Let q = z = const and 6 — const be material surfaces. Then B 
becomes 

B = VzxV8 = kx\78 = -V x . (1.32) 

So far the velocity field has been left u free but if this is then chosen 
such that in K 2 

u = -V^tp with 8 = -(-A) 1/2 V> (1.33) 

then the equations for B with this u satisfies 

DB 

— =BVu. (1.34) 

These are the 2D-QG equations discussed by Constantin, Majda & 
Tabak (1994) who linked the formation of a singularity to the a presence 
of hyperbolic saddle (for the level sets). Cordoba (1998) then showed the 
absence of a singularity in the case of a simple hyperbolic saddle. We 
refer the reader to for related work connected to the formation of sharp 
fronts and their evolution: Ohkitani & Yamada (1997) Constantin, Nie 
& Schorghofer (1998), Cordoba (1998), Cordoba, Fcfferman & Rodrigo 
(2004) & Rodrigo (2004). 



1.4 Transport equations for the curl and divergence of the 

Lamb vector 

For both the incompressible and compressible Euler equations, apart 
from the non-local effects of the pressure, the nonlincarity is the Lamb 
vector 

D = uxu. (1.35) 

This is the cross product of vorticity and velocity and is therefore an 
indicator of regions of a flow field where vorticity is nonzero. It shares 
its critical points with velocity and vorticity, but possesses additional 
critical points where the flow is locally Beltrami, as in the helical mo- 
tion of a swirling jet. The evolution of D is of interest : its divergence, 
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for example, plays a role in the production of jet-noise in compress- 
ible flows, wherever the mean of div D ^ 0. For recent discussions of 
the utility of the this vector as a diagnostic in fluid dynamics or as 
an important source of jet noise, see Rousseaux, Seifer, Steinberg & 
Wiebel (2007), Hamman, Klewick & Kirby (2008) and Cabana, Fortune 
& Jordan (2008), respectively. In the incompressible case, existence of 
solutions is assured provided the Beale, Kato, Majda (1984) condition 
is fulfilled. 

1.4-1 The evolution of D for the incompressible Euler 

equations 

In the case of the incompressible Euler equations, we distinguish the 
Lamb vector D — us x u from the Bernoulli vector 

E = D + V(p + \u 2 ) . (1.36) 

The incompressible Euler fluid equations (|1.3|) are connected to D and 
Eby 

d t u = -D-V{p+\u 2 ) = -E, (1.37) 

which vanishes for steady flows to create Lamb surfaces, reviewed, e.g., 
in Sposito (1997). The Bernoulli vector E is distinguished from the 
Lamb vector D by its divergence, in that 

div E = , while in general div D ^ , (1.38) 

although they both share the same curl, i.e., curli? = curl£) = vs. We 
now choose to rewrite the vorticity equation (|1.5p as 

3tu> + curl£: = 0, (1.39) 

and thereby remove any gauge freedom in the curl -1 operation. That 
is, we choose the Bernoulli gauge, in which curl _1 tcr = E. 

The aim of this section is to show that the curl of the Lamb vector 
vj = curl E plays a role similar to that of 13 in Theorem 11.2.11 and 
its divergence div D obeys a conservation equation that introduces an 
augmented transport velocity field, in the same spirit as in Haynes & 
Mclntyre (1987). 

The Euler fluid equations imply the following equation for the evolu- 
tion of the Lamb vector, 

d t D-uxzu = Exu>, (1.40) 
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and so 

Theorem 1.4.1 zu = curlE satisfies the stretching equation 

d t zz> - curl (uxtzt) = T>i am , (1.41) 
where T>i am is defined by 

T>i am = cm\{E x w), 
and div D satisfies the conservation equation 

d t (divD) +div[W(divD)] = 0, (1.42) 
where the transport velocity 11 is defined by 

(W-u)divD = «x(2S^)+ux V(p+ ±u 2 ) , (1.43) 
cmc? 5 is i/ie strain-rate tensor. 

Remarks : This theorem is similar in spirit to Theorem 11.2.11 for the 
evolution of 3 and the continuity equation for q, with za and div D 
playing these roles respectively. Four further observations about these 
equations follow that all hinge of the property that div E = 0. 

(1) Another expression for the divergence of the Lamb vector is 

div£> = - A{p+ \u 2 ). (1.44) 

Therefore, equation (|1.42[) is conceivably interesting as an evolution 
equation for the Bernoulli function (p+ \u 2 }. In compressible turbu- 
lence, such as in the exhaust of a jet airplane, the jet noise is largely due 
to correlations that produce a mean div_D ^ 0, as discussed in Cabana, 
Fortune & Jordan (2008). That is, the divergence of the Lamb vector 
is the leading source of turbulent jet noise, so the conservation equation 
(|1.42[) for its evolution in the incompressible case may be of interest. 
The quantity in square brackets in (|1.42[) is the current density for the 
transport of the hydrodynamic charge density, div_D = — A [p + \u 2 }. 

(2) Because &xv E = 0, the Helmholtz equation (|1.39[) and the curl of 
(|1.41[) . rewritten as 

zu t — curl (ti x zu) = curl (E X u>) , (1-45) 

imply the following two-component system of commutator equations 

zut + [u, zu] = [lj, curl - tp], (1-46) 
dtu + [u,b)] = 0, where u = curl -1 u; (1-47) 
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and we have E = curl -1 'n7 in the Bernoulli gauge. The bracket [• , •] 
in these equations denotes commutator of divergence-free vector fields. 
For example, 

[w, E] := uj -VE - E ■ Vw. (1.48) 

(3) Given that div E = 0, one may compute the evolution of E-helicity, 
defined as 



A B := 



J E-dxAd(E-dx) = J E - curl E dV . (1.49) 

Equations (|1.41[) and (|1.40[) imply that the helicity of E is not constant. 
Instead, evolves as 

d 



dt 



J E ■ curl EdV = -2 J u-(Ex curl E) dV , (1.50) 



after integrating by parts to remove gradient terms and applying homo- 
geneous boundary conditions. 

(4) Finally, we remark that the steps taken to distinguish between the 
fields D and E in selecting the Bernoulli gauge etc. are all reminiscent 
of a formal Eulerian analogy with Maxwell's equations for electromag- 
netism. All that remains in completing that well-known formal analogy 
is to identify u> with the magnetic field B and require the curl of the 
magnetic induction H to vanish, i.e., curlff = 0. Then one may inter- 
pret dtD in equation (|1.40[) as the displacement current, the right hand 
side as the current density, etc. . This is slightly different from the vari- 
ant of that formal analogy discussed previously in Marmanis (1998) and 
pursued in fluid experiments by Rousseaux, Seifer, Steinberg & Wiebel 
(2007). For completeness, we list the comparison in Table [T~2l 



Maxwell's equations 


Marmanis (1998) 


Present paper 


Magnetic Field, B 


w = curl u 


uj = curl u 


Magnetic Induction, H 


Absent 


Vx 


Electric Field, E 


LJ x u 


UJ X U + V(p+ ±M 2 ) 


Displacement vector, D 


Absent 


UJ X u 


Charge density, qE 


div (uj x u) 


div (uj x u) 



Fig. 1.2. Compared with Marmanis (1998), the current Maxwell- 
hydrodynamics analogy distinguishes between (E,B) and (D,H). 
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We shall refrain, however, from following the formal analogy between 
hydrodynamics and Maxwell equations here, and finish by interpreting 
the transport theorems for the Lamb vector's divergence (|1.42l) and its 
curl (| 1 .45[) in the standard fluid context. 

In the fluid context, we interpret the commutator equations (jl.46l) 
and (|1.47[) as evolution equations for the Lagrangian fluxes vj ■ dS and 
u> ■ dS as they are swept along by the fluid velocity u. Namely, 

j(vj-dS) = curl (E x w) • dS (1.51) 

|(wd5) = 0, (1.52) 

both along dx(t)/t = u(x(t),t). As always, the Helmholtz equation 
(11.52[) states that the flux of vorticity is frozen into the flow. However, 
the flux of the cross product of vorticity and the Bernoulli vector drives 
the flux of the the Lamb vector's curl, which in turn drives the evolution 
of the vorticity. In a nonlinear feedback response, the Lamb vector's curl 
is driven itself in equation (|1.5ip by the flux of the vector product of the 
vorticity with the Bernoulli vector, which contains both the nonlinearity 
and the pressure gradient. This process is akin to a magnetic dynamo, 
with the vorticity playing the role of the magnetic field. 

One may also express this process as a pair of linked circulation the- 
orems, namely, 

d 
dt 
d 

dt Jc(u) 

where c(u) is a closed material loop moving with the fluid velocity u. 
Thus, the circulation of the Lamb vector (or the Bernoulli vector) is 
driven by the circulation of the cross product of the Bernoulli vector 
with the vorticity. 

In addition, the Lamb vector's divergence (div D) satisfies the conti- 
nuity equation (| 1 .43(1 with its augmented transport velocity that depends 
on the vorticity, the strain-rate tensor and the pressure gradient. 



E ■ dx = <b (E x us) ■ dx (1.53) 

J c(u) 

dx = 0, (1.54) 



1-4-2 Helicity density 

Having looked at the dynamics of the Lamb vector D = u> x u let us 
consider the helicity density of the Euler equations which is the scalar 
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product 

A = uj ■ u . (1.55) 
Straightforward differentiation gives its dynamical equation, 

^- = -u,-V(p-\u 2 ) (1.56) 
which may be rewritten equivalently as 

<9 t A + div {Xu + uj (p- ±u 2 )} = 0. (1.57) 
As in tjl.21 this leads to the definition of a transport velocity field tlx 

\{U x -u)=u}{p-\u 2 ) (1.58) 
and the continuity equation 

d t X + dw{\Ux) = 0. (1.59) 
Thus, the vector quantity 

B x = VA x V0 (1.60) 
satisfies the stretching and folding result of Theorem 11.2.11 of £11.21 

d t &x-curl(ux&x) = Vx, (1.61) 
with vector T>x defined as 

V\ = -V(AdivW A ) x V0. (1.62) 

The vector T>x measures the "permeability" or rate of slippage of level 
sets of helicity density through level sets of the passive scalar field, 9. 

1.5 Conclusion 

The stretching and folding processes that produce small-scale structures 
in either fluid turbulence or MHD have generally been associated with 
the alignment or anti-alignment of either the vorticity u) or the magnetic 
field B with eigenvectors of the velocity gradient matrix Vit. The obser- 
vations and calculations in this paper have shown that these stretching 
and folding processes occur quite widely : on the one hand they have 
been shown to apply to any system that involves two passive scalars 
riding on a flow u, such as (q, 9) for the stratified, rotating Eulcr and 
Navier-Stokes equations, while on the other hand both the curl of the 
Lamb vector and the helicity density A = u> ■ u for incompressible flow 
also possess this behaviour. The significant feature is that we are one 



16 



gradient higher on u> itself. The embedded gradient of u) within 3 has 
been identified as the basis of a test for Euler codes, as explained in 

EB 

That the Lamb vector fits into this stretching and folding picture is 
a surprise. It is generally associated with studies in jet-noise in aero- 
acoustics and its natural context is compressible flows in which wave 
motion is observed, However, because it is the kinematic nonlinearity of 
fluid flow its evolution is important in also characterizing incompressible 
fluid motion. The stretching and folding process turns out to fit into 
an electro-magnetic analogy. We hope this analogy may become useful 
in transferring methods of mimetic difference schemes that are highly 
developed for electro-magnetic applications into the arena of Eulerian 
fluid dynamics. Mimetic methods are reviewed, for example, in Lipnikov, 
Shashkov & Yotov (2009). 

Acknowledgements: The authors thank Colin Cotter, Charles Doer- 
ing, Edriss Titi and Claude Bardos for helpful discussions. 
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